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I.  INTRODUCTION 


During  the  past  few  years,  there  have  been  many  advancements  In  the 
numerical  techniques  for  predicting  the  behavior  of  complex  fluid  flows. 
For  example,  several  computer  models  have  been  developed  by  Gosman, 
Spalding  and  others  [1,2,3]  which  use  the  mass,  momentum  and  energy 
conservation  equations  reduced  to  finite  difference,  nonlinear  algebraic 
form.  The  development  of  reliable  computer  programs  of  this  type 
greatly  benefits  engineering  analysis  in  such  widely  varying  fields  as 
meteorology,  aerodynamics  and  gasdynamics. 

The  earlier  two-dimensional  computer  codes  were  based  on  vorticity 
(w)  and  stream  function  (<|/)  [1,2,4].  This  form  of  the  governing  equa¬ 
tions  eliminates  pressure  and  velocity  from  immediate  consideration. 
Pressure  is  normally  calculated  only  after  a  converged  solution  is 
obtained.  This  technique  has  several  inherent  disadvantages: 

1.  It  results  in  large  errors  in  the  predicted  pressure  distribu¬ 
tions  in  all  but  quiescent  flow  regions  due  to  the  higher  order  depen¬ 
dence  of  the  pressure  gradient  on  stream  function  [5]. 

2.  It  is  usually  restricted  to  constant  density  flows  or  to  flows 
in  which  density  varies  only  with  temperature  [3,5]. 

3.  The  boundary  conditions  are  difficult  to  specify  [3,5]. 

4.  Considerable  difficulty  has  been  experienced  in  arriving  at 
converged  solutions,  especially  for  nonuniformiy  spaced  grids  and  high 
flow  rates  [2,4,5] . 

5.  The  model  is  not  easily  extended  to  three  dimensional 
flows  [3]  • 
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To  overcome  these  difficulties,  emphasis  has  been  placed  on  developing 
computer  codes  based  on  velocity  and  pressure,  the  primitive  variables. 

A  major  problem  with  any  new  computer  model  is  model  validation. 

The  difficulties  of  collecting  accurate  empirical  data  are  multiplied 
wherv  investigating  three  dimensional  and/or  reacting  flows.  In  addition, 
many  variables  within  these  flows  are  not  readily  measurable  (turbulence 
intensities,  etc.). 

An  effort  to  utilize  elliptic  computer  models  which  can  handle 
turbulent,  reacting,  variable  density  flows  at  high  subsonic  and  sonic 
velocities  has  been  underway  at  the  Naval  Postgraduate  School  for  several 
years.  Two  specific  areas  which  have  been  Investigated  are  flows  in  a 
turbojet  test  cell  and  in  the  combustion  environment  of  a  solid  fuel 
ramjet. 

A  solid  fuel  ramjet  (SFRJ)  most  often  consists  of  a  solid  fuel  grain 
which  provides  the  walls  for  the  combustion  chamber  [A].  Located  at  the 
air  inlet  end  of  the  combustor  is  a  sudden  expansion  or  other  type  of 
flame  stabilization  device.  The  opposite  end,  downstream  of  the  fuel 
grain,  may  also  incorporate  a  sudden  expansion  aft  mixing  chamber.  The 
primary  combustion  region  contains  a  turbulent  diffusion  flame  which  emanates 
from  the  forward  recirculation  zone  and  remains  within  the  developing 
boundary  layer.  The  aft  mixing  region  may  incorporate  some  means  of 
injecting  air  (bypass  air)  in  order  to  complete  the  consumption  of  the 
fuel  which  exits  the  aft  end  of  the  fuel  grain.  Mixing  chamber  and  inlet 
design  variables,  fuel  grain  design  and  fuel  properties  make  a  wide  variety 
of  performance  characteristics  available. 


2 


The  possibility  of  incorporating  this  type  of  propulsion  device  into  a 
future  medium  or  long-range  tactical  weapon  system  coupled  with  the  ex¬ 
pense  of  testing  each  new  design,  makes  the  continued  development  of 
reliable  computer  models  highly  desirable.  The  model  could  be  used  to 
predict  the  effects  of  fuel  properties  and  to  inexpensively  evaluate 
different  geometries  and  operating  conditions.  In  addition,  a  three 
dimensional  code  would  allow  modeling  discrete  air  injection  into  the  aft 
mixing  region.  The  latter  technique  can  substantially  increase  combus¬ 
tion  efficiency  and  allowable  fuel  loading. 

Previous  work  at  the  Naval  Postgraduate  School  has  been  directed 
toward  improvement  of  the  quantitative  accuracy  of  the  model  and 
toward  validation  of  that  model  [A] .  Reasonable  agreement  with  empirical 
data  has  been  obtained.  However,  as  previously  stated,  the  ip-u  model 
does  not  predict  accurate  pressure  distributions  and  numerical  difficul¬ 
ties  prevented  modeling  the  aft  mixing  chamber. 

The  purpose  of  this  investigation  was  to  adapt  and  validate  a  primi¬ 
tive  variable,  two-dimensional,  finite  difference  computer  code  which 
models  the  flow  within  a  solid  fuel  ramjet. 
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II.  MODEL  OVERVIEW 


A.  INTRODUCTION 

The  computer  model  used  in  this  study  was  adapted  from  the  CHAMPION 
2/E/MX  computer  program  developed  by  Pun  and  Spalding  [6].  CHAMPION  is 
a  TVO~dimensional  Elliptic,  FIXed  grid  computer  program  which  provides 
a  solution  of  the  conservation  equations  for  recirculating  flows  in 
finite  difference  form. 

B.  ASSUMPTIONS 

The  flow  was  assumed  to  be  steady,  two-dimensional  and  subsonic. 

For  simplicity  the  value  of  specific  heat  (C^)  was  assumed  to  be  constant 
although  its  dependence  on  temperature  and/or  composition  could  easily 
be  included. 

A  modified  Jones-Launder  [6,7,8,91  two  parameter  turbulence  model 
was  incorporated  to  calculate  the  effective  viscosity.  It  uses  five 
empirical  constants  (Table  I)  and  requires  that  two  additional  variables, 
turbulencekinetic  energy  (K)  and  turbulencedissipation  rate  (e) ,  be 
evaluated.  Effective  viscosity  was  calculated  using  the  formulas: 

yeff  "  Wlam  +  yt  (,) 

where 

yt  -  CD  p  K2/e  (2) 
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For  reacting  flows,  the  four  species,  oxygen,  nitrogen,  fuel,  and 
products,  were  considered.  Simple,  one-step,  infinitely  fast  kinetics 
were  assumed  in  which  a  fuel  combines  with  an  oxidant  to  form  a  single 
product  without  intermediaries  [4,10]. 

1  kg  fuel  +  i  kg  oxidizer  -►  (1  +  i)  kg  products 

Fuel  and  oxygen,  therefore,  could  not  exist  simultaneoulsy  and  the  com¬ 
bustion  process  was  mixing  limited,  in  addition,  it  was  assumed  that  no 
oxygen  existed  at  the  fuel  surface  and  that  surface  was  isothermal.  The 
turbulent  Prandtl  and  Schmidt  numbers  were  taken  equal  to  unity  and, 
therefore,  the  turbulent  Lewis  number  was  unity.  The  laminar  Prandtl 
number  was  also  taken  to  be  unity. 

C.  GOVERNING  EQUATIONS 

The  conservation  equations  for  axi -symmetri cal  flows  with  no  tangential 
variations  can  be  put  into  the  general  form  [6]: 

|^pu(0)  -•--l^prv*)  -  vrf^riylf)  -  (3) 

convection  terms  diffusion  terms  source  terms 

where  d>  stands  for  the  dependent  variable  (u,v,k,£,h,  etc..)  being  con¬ 
sidered  ($  *  1  for  the  continuity  equation),  is  the  appropriate  effec¬ 
tive  exchange  coefficient  for  turbulent  flow  and  is  the  "source  term" 
(Table  ll).  The  energy  equation  in  terms  of  stagnation  enthalpy  has  no 
source  terms  since  the  turbulent  Prandtl  and  Schmidt  numbers  were  chosen 
as  unity  and  radiative  transport  was  neglected  [1,3].  Thus  the 
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stagnation  enthalpy  is  given  by: 


h  -  h  +  (u2  +  v2)/ 2  +  K 


(*) 


where  for  non-reacting  flows: 


h  =  C  T 
P 


and  for  reacting  flows: 


h  =  mox  AH/I  ♦  Cp(T  -  Tref) 


(5) 


(6) 


The  calculation  of  temperature  was  made  using  equations  (4),  (5) 
and  (6).  Density  was  calculated  from  the  perfect  gas  law: 

p  -  P/RT 

# 

for  non-reacting  flows:  R  *  constant 
for  reacting  flows:  R  *  R/M 


where, 


i  / •*  I,,  +  m  /M  ^  ^  m  /M 

1/M  =  m-  /M,  T  mox  ox  N2  N2  pr  pr 
'  u  ru 


(7) 


(8) 


(9) 


For  modeling  reacting  flows,  two  additional  quantities,  m^  and 
X  =  m^u  -  ,  were  evaluated.  Each  of  these  properties  as  well  as 

stagnation  enthalpy  have  identical  governing  differential  equations 
(equation  (3)  with  no  source  terms),  in  appropriate  dimensionless  form 
they  also  have  identical  boundary  conditions.  Thus,  only  one  of  the 
equations  had  to  be  solved.  The  dimensionless  form  selected  for  each 
property  was: 


H  -  <hin  '  «'<hm  '  hfg> 


(10) 
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l2)/("M2,n  '  %2fs> 

(ID 

■t,  - 

(12) 

In  this  study,  stagnation  anthalpy  was  calculatad.  H  was  then  formed 
using  equation  (10).  Since  H  ■  m^  ■  X  all  points  in  the  flow  field, 
1^2  and  x  could  be  calculated  using  equations  (11)  and  (12).  The  mass 
fractions  of  fuel,  oxygen,  and  products  (m^u,  f»ox,  mpr)  were  found  from 
the  equations: 

for  x  >  o;  mfu  *  X.  mox  ■  0 


for  x  <  0;  mfu  -  0,  mQx  -  -xi 


mpr  '  1  '  "|>r  -  V  '  mfu  '  "W  (,4) 

D.  CONSERVATION  OF  MASS 

On  each  radial  line  the  mass  flow  rate  was  calculated  using  the  local 
density.  The  error  in  mass  flow  (compared  to  the  summation  of  "mass-in" 
at  all  upstream  boundaries)  was  used  to  uniformly  adjust  the  axial  veloci¬ 
ty  over  the  entire  line.  This  process  ensured  that  overall  continuity  was 
satisfied  on  the  line.  The  pressure  at  all  downstream  locations  was  then 
adjusted  to  approximately  correct  for  the  momentum  imbalance  created  by 
the  uniform  axial  velocity  adjustment.  A  "pressure  correction"  equation 
was  then  solved  for  each  cell  on  the  line.  Local  cell  velocity  (axial 
and  radial)  and  pressure  were  then  adjusted  to  satisfy  cell-wise  conti¬ 
nuity.  The  details  of  this  procedure  are  presented  in  reference  6. 
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E.  BOUNDARY  CONDITIONS 


1 .  Introduction 

Fixed  boundary  conditions  were  specified  at  the  desired  or  ex¬ 
perimentally  determined  values.  Specified  gradient  boundary  conditions 
were  handled  by  setting  the  appropriate  convection/diffusion  coefficient 
to  zero  in  the  finite  difference  equation  ("breaking  the  link")  and  then 
entering  the  appropriate  gradient  through  linearized  "false"  source  terms 
[6].  The  geometry  and  appropriate  boundary  conditions  are  summarized  in 
figure  1. 

2.  Inlet 

Although  not  a  computer  program  limitation,  "plug  flow"  was  assumed 
at  the  inlet.  turbulence  kineti c  energy  was  selected  to  be  uniform  with  a 
value  which  corresponded  to  the  approximate  turbulence  intensity  of  the 
inlet  flow. 

3.  Axis  of  Symmetry  and  Exit  Plane 

Radial  and  axial  gradients  were  set  equal  to  zero  on  the  center 
line  and  exit  respectively.  The  radial  flow  velocity  was  equated  to  zero. 

4.  Solid  Boundaries 

All  non-reacting  solid  boundaries  were  considered  adiabatic  with 
both  velocity  components  equal  to  zero  ("no  slip"  condition). 

For  simplicity,  a  two  part  boundary  layer  was  used.  The  border 
between  the  laminar  sublayer  and  the  turbulent  layer  was  taken  at 

yp+  «  11,5  [6],  yp+  was  evaluated  at  each  near  wall  node  (p) , 

C*  1/2  (l5) 


(16) 


where,  for  y  >11.5 
P  — 


Tw  -  CD,/2  p  Kp 


Tw  was  assumed  uniform  from  the  wall  to  the  near  wall  grid  point.  Thus, 


♦  ■  p  K  1/2  S/ll 


»P 


1am 


(17) 


If  y  +  11.5,  the  wall  shear  stress  (t  )  was  calculated  using  the  formula 

p  1  w 


Tw  *  CD1/2  P  Kp  "  P  CD,/4  <p1/2(u/u+) 

“  *  Co’A  p  Up  Kp,/2/ln(EP6CD1/\,/2/Plam> 


(18) 


where 


u+  £  ,n^E  yp+^ 


(19) 


Wall  shear  stress  was  evaluated  for  yp  <11.5  from  the  formula: 


w 


“lam  Vs 


(20) 


Due  to  the  steep  gradients  of  properties  in  turbulent  flows  near 
solid  boundaries,  the  source  terms  for  K  and  e  at  near  wail  nodes  were 
expressed  in  terms  of  the  wall  shear  stress  [1,6].  also  provides  the 
boundary  condition  for  the  u  and  v  equations.  In  the  following  equation 
for  turbulence  dissipation  rate  (e)  at  a  near  wall  node  (p) ,  the  length 
seal*  is  presumed  proportional  to  the  distance  from  the  wall  (<5) . 

Ep  -  C03/A  Kp3/2/K«  -  Kp3/2/2.43<5  (21) 
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It  was  found,  as  was  previously  found  by  Netzer  [4],  that  (when  using  the 

sudden  expansion  geometry  in  reacting  flows)  the  near  wall  dissipation 

3/2 

had  to  be  increased  on  the  step  face  (e  ■  K  3  /0.k6)  and  that  the  grid 

P  P 

spacing  adjacent  to  the  fuel  surface  had  to  be  fine  (yp+  <  11.5)  in 
order  to  obtain  a  temperature  distribution  in  qualitative  agreement  with 
experiment.  £quation  (20)  implies  that  the  wall  shear  stress  is  calcula¬ 
ted  assuming  a  linear  velocity  profile  when  yp+  <11.5.  A  near-wall 
grid  point,  therefore,  can  lie  within  the  laminar  sublayer,  but  the 
source  terms  for  K  and  e  imply  that  Is  much  greater  than  one 

[7,81.  This  fact  precludes  yp+  from  being  significantly  less  than  11.5. 

For  reacting  flaws,  the  boundary  conditions  for  the  dimensionless 
properties  (equations  (10),  (li)  and  (12))  were  zero  at  the  inlet  and 
unity  "deep"  in  the  fuel  grain  (fg).  These  properties  were  considered 
to  have  zero  gradients  on  non- reacting  surfaces. 

The  assumptions  employed  for  reacting  flows  (unity  turbulent  Prandtl 
and  Schmidt  numbers,  simple  chemical  reaction,  constant  specific  heat  and 
stagnation  enthalpy  defined  in  equations  (4)  and  (6))  result  in  a  general 
boundary  condition  for  all  "conserved"  properties  (<J>c)[10]  on  a  surface 
which  has  mass  transfer, 

3<fc 

*S-  *  sAw'V  -  V 

where  <b  represents  h,  mu,  or  y  =  m*,,  “in  /i  • 

C  Na  t  u  ox 

A  mass  transfer  cnnductance  (g)  is  often  defined  such  that, 

3<t> 

<r*  3T*bw  ’  123 

where  $  is  defined  as  the  free  stream  value.  For  this  application,  d> 


was  taken  to  be  the  local  near  wall  value  4  . 

P 

Substituting  equation  (23)  into  equation  (22)  yields: 


(24) 


"’"bw  ■  "  *c  >/<*e  ■  * 

***  cp  Cbw  Cbw  Cfg 


=  g  BP 


(25) 


where  BP  represents  the  mass  transfer  (or  "blowing")  parameter. 

Without  mass  transfer  the  wail  heat  flux  (qw")  can  be  defined  in 
terms  of  the  conditions  at  the  near  wail  node. 


V  <»> 

where  hw  is  the  enthalpy  of  the  wall  and  h  is  the  heat  transfer  conduc¬ 
tance. 

With  <p  »  h  in  equation  (23), 
c 

9’(^!?')„/(hp-hw)  (27> 

Substi tuting  equation  (26)  into  equation  (27)  yields: 

9  -  h/ Cp  (28) 

or 

g/(pu)p  -  fc/[(pu)p  cp]  =  St  (29) 

thus, 

g  ■  (pu)p  st  (30) 

From  Reynolds  Analogy  with  unity  Prandti  number, 

St  -  Cf/2  -  tw/(pu2)p  (31) 

where  Cf  is  the  local  friction  coefficient. 
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Combining  equations  (30)  and  (31)  yields: 

9  -  Tw/up  (32) 

Using  the  Couette  flow  approximation  for  the  boundary  layer  behavior 
with  mass  transfer  [10], 

g  -  g*  1 n ( 1  +  BP)/BP  (33) 

where 

g*  '  ,5mBP-K>^  (31*) 

In  this  application,  BP  was  evaluated  from  the  solution  of  the 
energy  equation  using. 


BP  -  (hp  -  hJ'V  -  hf9)  (35) 

The  wall  shear  stress  was  calculated  using  equation  (18)  or  equation 
(20)  and  modified  with  equation  (33). 

Tbw"  Tw  ,n0  +  BP)/BP  (36) 

where  t  is  the  wall  shear  stress  without  wall  mass  addition, 
w 

The  mass  transfer  conductance  (g)  was  found  using  equation  (32).  The 
wall  mass  flux  was  then  evaluated  using  equation  (25). 

The  wall  heat  flux  (dw")  on  all  solid  isothermal  boundaries  was  evalu¬ 
ated  using  the  Reynolds  analogy: 

-V'^p  -  \,>  -  V“P  (37) 
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Since  the  blowing  rates  were  smell  for  the  solid  fuel  ramjet  (typi¬ 
cally,  BP  <  2.0),  K  and  e  were  evaluated  using  equation  (3)  and  the  terms 
presented  in  Table  II  which  incorporate  the  empirical  constants  of  Table 
I. 

Blowing  velocity  (v^w)  and  fuel  regression  rate  (RR)  were  calcula¬ 
ted  using  the  formulas: 

vb*  *  •|4bw"/(Ibw  <3®> 

RR  ■  (39> 

F.  SOLUTION  PROCEDURE 

Five  variables  (u,v,K,e  and  H  or  h)  were  solved  using  equation  (3)  in 
finite  difference  form.  The  line  by  line  iterative  procedure  employed 
upwind  differencing  and  under  relaxation  to  promote  convergence  [6]. 
Pressure  (relative  to  a  selectable  position  and  magnitude  within  the 
grid)  was  obtained  from  the  mass  conservation  imposed  on  each  radial 
grid  line  and  on  each  nodal  control  volume  as  discussed  above.  Effective 
viscosity,  temperature  and  density  were  also  obtained  as  described  above. 

A  more  detailed  explanation  of  this  procedure  can  be  found  in  reference6. 


III.  DISCUSSION  OF  RESULTS 


A.  INTRODUCTION 

The  purpose  of  this  study  was  to  develop  a  primitive  variable,  finite- 
difference  computer  program  that  could  be  used  to  determine  the  flow  within 
a  solid  fuel  ramjet  combustor  wl th  emphasis  on  the  aft  mixing  chamber.  The 
effects  of  inlet  mass  flow  rate  and  inlet  dump  area  ratio  on  the  flow  field 
were  examined.  As  previously  explained,  an  aft  mixing  region  allows  further 
combustion  aft  of  the  fuel  grain.  This  process  normally  increases  combustion 
efficiency.  Lowering  the  inlet  flow  rate  increases  the  fuel-air  ratio  within 
the  fuel  port.  Bypass  air  can  then  be  injected  into  the  aft  mixing  region. 

The  latter  procedure  can  be  used  to  appreciably  increase  fuel  loading.  Pre¬ 
vious  work  at  the  Nava)  Postgraduate  School  [4,11,12]  modeled  a  SFRJ  with  a 
computer  program  utilizing  <J/-u  as  primary  variables.  Numerical  instabilities, 
however,  prevented  the  use  of  the  ty-w  model  to  predict  the  flow  in  the  aft 
mixing  region.  The  results  of  that  investigation  and  some  empirical  data  were 
available  for  comparison  with  the  predictions  from  the  primitive  variable 
model. 

Several  factors  were  anticipated  which  could  contribute  to  differences 
between  the  predictions  of  the  two  models  and  the  empirical  data: 

a.  Some  of  the  experimental  data  were  measured  in  cold,  nonreacting 
flows. 

b.  The  incorporation  of  the  aft  mixing  chamber  into  the  primitive 
variable  model  could  influence  the  flow  upstream  in  the  combustion  chamber. 

c.  In  the  model,  a  wall  value  of  turbulence  kinetic  energy  (K) 

was  specified  through  a  slip  factor  such  that  K  ■  (-1.0  or  -0.39)*K_, 

w  p 
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depending  on  the  magnitude  of  the  turbulent  Reynolds  Number.  In  the 
u-v-p  model  the  boundary  condition  for  K  at  the  near  wall  node  (p)  was 
specified  in  terms  of  the  wal 1  shear  stress.  In  addition,  in  the  primi¬ 
tive  variable  model,  the  boundary  condition  for  stagnation  enthalpy  at 
the  near  wall  node  was  made  a  function  of  wall  shear  stress  through  the 
Reynolds  Analogy.  These  factors  affect  heat  flux  to  the  wal l ,  and, 

therefore,  the  fuel  regression  rate. 

d.  The  u-v-p  model  incorporated  a  23  by  21  grid  in  the  fuel 
port  while  the  model  utilized  a  17  by  25  grid.  In  reality  the 
heat  of  vaporization  of  the  fuel  is  a  fixed  quantity  and,  if  converged 
solutions  are  obtained,  the  wall  heat  flux  should  not  depend  upon  the 
grid  spacing.  However,  it  has  been  found  [4]  that  the  heat  flux  to  the 

wall  (which  is  calculated  using  the  near-wall  grid  point)  is  a 

function  of  the  grid  distance  from  the  wall.  This  results  from  the 
assumed  behavior  of  the  variables  near  the  wall.  The  procedure  employed 
in  this  study  was  to  adjust  the  heat  of  vaporization  to  match  the  empiri¬ 
cal  fuel  regression  rate  at  one  air  flow  rate  and  to  use  that  value  for 
all  other  flow  rates.  If  the  model  is  realistic,  fuel  regression  rate 

should  then  vary  with  air  flow  rate  in  agreement  with  experiment. 

8.  REGRESSION  RATE 

Figure  2  shows  that  the  fuel  regression  rate  predictions  of  the  u-v-p 
and  ip-oi  models  are  quite  similar.  Both  predict  the  peak  regression  rate 
upstream  of  experiment  and  have  similar  slopes.  This  early  peak  in  the 
regression  rate  results  from  the  model  predicting  a  shorter  reattachment 
length  than  was  found  experimentally  [4].  The  primitive  variable  model 
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predicted  higher  regression  rates  downstream  of  flow  reattachment  in 
better  agreement  with  experiment. 

Figure  3  shows  the  effects  of  increasing  inlet  air  mass  flux  (G  ■ 
majr/Aj.p  a  Vjn)  on  fuel  regression  rate  (r^u).  The  regression  profile 
remained  the  same  and,  as  expected,  decreased  with  decreasing  G.  It 
has  been  found  experimentally  that  the  regression  rate  of  plexiglass 
varies  as  the  air  mass  flux  raised  to  a  constant  power  (r^  <*  Gn) .  Boaz 
and  Netzer  [11]  found  that  this  constant  was  equal  to  0.4l  while  Mady, 
et  al  [12]  found  it  to  be  approximately  0.38.  For  the  three  test  cases 
of  this  study,  the  u-v-p  model  predicted  the  constant,  n,  to  be  between 
0.31  and  0.34.  Thus,  the  primitive  variable  model  appears  to  correctly 
predict  the  nature  of  the  change  in  convective  heat  flux  to  the  fuel 
surface  with  air  flow  rate. 

C.  TURBULENCE  INTENSITY 

Figure  4  comparies  the  predicted  centerline  turbulence  intensity 
(assuming  isotropic  turbulence)  and  experimental  data  for  non-reacting 
flow.  The  primitive  variable  computer  model  slightly  underpredicted  the 
peak  turbulence  intensity  while  the  iJj-u)  model  overpredicted  it.  Both 
models  predicted  the  peak  occurring  downstream  of  experiment  and  both 
distributions  appear  to  approach  an  identical  asymptote  downstream. 

The  decrease  in  turbulence  intensity  predicted  by  the  ty-w  model  near  the 
inlet  resulted  from  the  model  over-predicting  the  velocity  increase  as 
the  air  entered  the  combustor  [4].  The  u-v-p  model  overcame  this  diffi¬ 
culty.  The  differences  in  the  results  from  the  two  computer  models  may 
result  from  the  differences  in  the  boundary  conditions  on  turbulence 
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kinetic  energy  in  the  combustor  and/or  to  the  effects  of  the  addition 
of  the  aft  mixing  chamber  on  the  upstream  flow.  It  should  also  be  noted 
that  the  experimental  data  used  in  this  comparison  were  obtained  in  a 
non- reacting  flow. 

Figure  5  shows  the  effect  of  decreasing  inlet  air  mass  flux  on 
turbulence  intensity.  As  anticipated,  the  peak  turbulence  intensity 
decreased  as  inlet  axial  velocity  decreased.  Each  test  condition,  however, 
converged  on  the  same  value  downstream.  Much  additional  experimental  work 
is  required  to  obtain  the  turbulence  intensities  in  reacting  flows;  only 
then  can  the  adequacy  of  the  K-e  turbulence  model  be  fully  evaluated. 

0.  PRESSURE 

Figure  6  shows  the  effect  of  inlet  velocity  on  the  axial  pressure 
distributions  for  the  three  primitive  variable  test  condi tions  (Table  III). 

(The  radial  scale  has  been  expanded  to  illustrate  the  pressure  variations. 

The  maximum  pressure  variation  is  approximately  1.2  psi.)  The  radial 
location  of  these  distributions  is  given  as  a  fraction  of  the  fuel  port 
radius  (R^p).  As  expected,  pressure  initially  increased  due  to  jet 
spreading.  This  was  followed  by  a  slight  pressure  drop  as  the  flow  ac¬ 
celerated  due  to  heat  addition  and  wall  friction.  The  final  pressure 
rise  was  due  to  jet  spreading  In  the  aft  mixing  region. 

E.  TEMPERATURE 

Figure  7  displays  radial  temperature  variations  in  the  combustor 
near  the  end  of  the  full  grain  and  at  about  1.5  aft  mixing  region  dia¬ 
meters  down  the  aft  chamber.  As  discussed  above,  fuel  flow  rate 
decreases  less  than  inlet  air  flow  rate  (rf  “  Gn,  n  <  1).  Therefore, 


as  air  flow  rate  is  decreased,  the  overall  mixture  ratio  becomes  more 
fuel  rich,  and  the  developing  boundary  layer  and  the  fuel  layer  between 
the  diffusion  flame  and  the  wall  thicken  .  Thus,  as  shown  in  figure  7,  as 
the  inlet  velocity  (and,  therefore,  the  inlet  air  mass  flux)  was  de¬ 
creased,  the  maximum  temperature  (or  "flame11)  in  the  combustor  moved  away 
from  the  fuel  grain  and  the  centerline  temperature  increased.  The 
maximum  temperature  in  the  aft  mixing  chamber  was  also  predicted  to  occur 
farther  from  the  top  wall. 

Figure  8  shows  similar  data  predicted  by  the  model  slightly 
farther  upstream.  A  significant  difference  between  the  predictions  of 
the  two  computer  models  was  that  the  model  predicted  a  stronger 
dependence  of  the  peak  temperature  radial  location  on  the  inlet  air  ve¬ 
locity.  An  aft  mixing  region  was  not  incorporated  into  the  ip-m  model. 
Therefore,  the  boundary  layer  continued  to  grow  and  the  point  of  maximum 
temperature  continued  to  recede  from  the  fuel  surface  with  increasing 
axial  distance  from  the  initial  reattachment  point.  The  aft  mixing  region 
of  the  u-v-p  model  caused  the  boundary  layer  thickness  (and,  therefore, 
the  location  of  the  peak  temperature)  to  become  approximately  constant  in 
the  latter  portion  of  the  combustion  chamber.  This  was  the  apparent 
cause  of  the  weaker  dependence  predicted  by  the  u-v-p  model  of  peak 
temperature  location  and  boundary  layer  thickness  on  inlet  air  mass  flux. 

Figure  9  is  an  illustration  of  the  predicted  combustion  behavior  in 
the  aft  mixing  region.  (The  radial  dimension  has  been  expanded  for  clarity.) 
Lines  of  maximum  temperature  (i.e.,  the  flame  sheet  location)  are  pre¬ 
sented  as  a  function  of  fuel  grain  inlet  air  velocity.  It  should  be  noted 
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that  the  aft  reel rculation  zone,  which  is  aiso  depicted  on  this  figure, 
was  predicted  to  be  fuel  rich  and  did  not  vary  appreciably  in  size  with 
changing  inlet  air  mass  flux.  As  discussed  above,  the  fuel  regression 
rate  decreased  more  slowly  than  inlet  air  flow  rate.  Thus,  as  air 
flux  was  decreased  the  mixture  entering  the  aft  chamber  became  more  fuel 
rich  and  the  thickness  of  the  fuel  layer  at  the  end  of  the  fuel  grain 
increased  slightly.  With  high  air  mass  flux  through  the  fuel  port  the 
mixture  ratio  is  fuel  lean.  The  flame  therefore  propagates  to  the  outer 
wall  of  the  aft  mixing  chamber.  This  condition  could  be  expected  to 
produce  a  high  combustion  efficiency,  it  should  be  noted  that  an  ade¬ 
quate  length-to-diameter  ratio  is  required  to  allow  the  flame  to  spread 
to  the  wall.  This  ratio  is  apparently  a  function  of  the  fuel  port  to  aft 
mixing  chamber  cross-sectional  area  ratio-  As  the  air  flow  rate  was 
decreased  the  mixture  ratio  became  fuel  rich  and  the  flame  did  not  reach 
the  wall*  This  would  result  in  unburned  fuel  entering  the  nozzle  and  a 
lower  combustion  efficiency.  Figure  10  shows  the  predicted  effect  of  the 
fuel  inlet  dump  step  size  on  the  flame  behavior  in  the  aft  mixing  cham¬ 
ber.  The  recirculation  zone  did  not  change  in  size  since  the  air  mass 
fluxes  were  identical.  The  smaller  inlet  step  produced  a  slightly  higher 
fuel  regression  rate  and  therefore  required  a  longer  aft  mixing  region. 
These  predictions  might  be  used  as  a  first  approximation  for  predicting 
the  "best"  placement  of  bypass  air  dumps  in  the  aft  mixing  region.  To 
predict  an  optimum  location,  however,  the  primitive  variable  model  would 
have  to  be  expanded  to  three  dimensions. 
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F.  COMPUTER  RELATED  PROBLEMS 


As  has  been  discussed  previous ly,  in  order  to  obtain  resuits  that  were 
in  agreement  with  experiment,  the  grid  spacing  near  the  fuel  surface  was 
required  to  be  fine  and  the  length  scale  of  turbulence  was  decreased  on 
the  combustor  step  face.  Because  convergence  was  sensitive  to  the  length- 
to-width  ratio  of  individual  control  volumes,  the  small  radial  grid  spacing 
near  the  fuel  surface  forced  similar  fine  spacing  in  the  axial  direction 
downstream  in  the  aft  mixing  region.  A  length  to  width  ratio  of  less  than  ten 
to  one  was  required.  These  criteria  forced  the  use  of  a  large  number  of  cells 
which  in  turn  required  a  large  amount  of  CPU  time.  A  typical  primitive 
variable  40  by  33  grid  required  75  to  80  minutes- of  CPU  time  on  an  IBM  360-67 
computer  to  converge.  A  typical  model  with  a  17  by  25  grid  required  35 
to  40  minutes  of  CPU  time.  It  must  be  remembered,  however,  that  numerical 
instabilities  prohibited  the  modeling  of  an  aft  mixing  chamber  with  the  ip-u> 
model . 

The  primitive  variable  model  demonstrated  some  convergence  difficulty  in 
the  aft  recirculation  region.  This  problem  seemed  to  be  associated  with 
the  continually  changing  velocity  profile  just  prior  to  the  aft  expansion 
(the  "inlet"  conditions  for  the  aft  mixing  chamber).  This  effect  was 
suppressed  by  sweeping  through  the  entire  flow  field  several 
times  with  only  a  few  traverses  on  each  line  and  then  increasing  the  number 
of  traverses  on  the  radial  grid  lines  in  the  aft  mixing  region  once  the 
combustor  flow  field  had  essentially  converged. 


IV.  CONCLUSIONS 


In  general,  the  predicted  flow  fields  for  the  two  computer  models 
were  quite  similar  within  the  fuel  grain.  However,  the  presence  of  the 
aft  mixing  region  coupled  with  the  few  boundary  condition  differences 
previously  mentioned,  had  some  effect  on  the  flow  field  predictions.  Th< 
most  noticeable  of  these  was  the  decrease  in  dependence  of  the  boundary 
layer  thickness  and  the  maximum  temperature  radial  location  on  axial  in¬ 
let  velocity.  As-anticipated,  the  primitive  variable  model  allowed  the 
prediction  of  the  flow  field  within  the  aft  mixing  region.  This  was  not 
possible  with  the  ip~w  model.  Many  additional  empirical  data  are  needed 
to  completely  assess  the  validity  of  the  primitive  variable  model  in 
predicting  the  flow  in  a  SFRJ. 
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FIGURE  3.  EFFECTS  OF  AIR  MASS  FLUX  ON  THE  PREDICTED  REGRESSION 
RATES  OF  THE  PRIMITIVE  VARIABLE  COMPUTER  MODEL 
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FIGURE  5.  PREDICTED  CENTERLINE  TURBULENCE  INTENSITY  AS  A  FUNCTION 
OF  INLET  VELOCITY  (PRIMITIVE  VARIABLE  COMPUTER  MODEL, 


FIGURE  6 .  PREDICTED  AXIAL  PRESSURE  DISTRIBUTION  AS  A  FUNCTION  OF 
INLET  VELOCITY  (PRIMITIVE  VARIABLE  COMPUTER  MODEL) 


FIGURE  7  .  PREDICTED  COMBUSTOR  AND  AFT  MIXING  CHAMBER  RADIAL  TEMPERATURE 
VARIATIONS  (PRIMITIVE  VARIABLE  COMPUTER  MODEL) 


FUEL  SURFACE 


FIGURE  o  .  PREDICTED  RADIAL  TEMPERATURE  DISTRIBUTIONS 

( ip— a)  COMPUTER  MODEL,  FIGURE  8  OF  REFERENCE  4) 


FIGURE  9.  Effect  of  Air  Flow  Rate  on  Predicted  Flame  Location  and 
Flow  Reattachment  Point  in  the  Aft  Mixing  Region  • 


the  Aft  Mixing  Region 


V.  REFERENCES 


1.  Gosman,  A.O. ,  and  others.  Heat  and  Mass  Transfer  in  Recirculating  Flows, 
Academic  Press,  1969. 

2.  Hayes,  J.D.  and  Netzer,  D.W. ,  An  investigation  of  the  Flow  In  Turbojet 
Test  Cells  and  Augmenters,  Naval  Postgraduate  School  Report  Number  NPS- 
57Nt-ftiO I,  October  1975. 

3.  Mongla,  H.C.  and  Reynolds,  R.S..  Combustor  Design  Criteria  Validation 
Volume  111,  Al RESEARCH  Manufacturing  Company  of  Arizona  Report  Number 
USARTL-TR-78-55C,  February  1979. 

4.  Netzer,  D.W. ."Model ing  Solid-Fuel  Ramjet  Combustion," Journal  of  Space¬ 
craft  and  Rockets,  Volume  14,"  Number  12,  p.  762-766,  December  1977. 

5.  Walters,  J.J.  and  Netzer,  D.W. ,  A  Validation  of  Mathematical  Models  of 
Turbojet  Test  Cells,  Naval  Postgraduate  School  Report  Number  NPS67-78-002, 
June  i 978 . 

6.  Pun,  W.M.  and  Spalding,  D.B.,  A  General  Computer  Program  for  Two-Dimensional 
El  1 iptic  Flows,  Imperical  College  of  Science  and  Technology,  Report  No. 
HTS/7672,  August  1977. 

7.  Launder,  B.E.  and  Spalding,  D.B.,"The  Numerical  Computation  of  Turbulent 
Flowd;  Computer  Methods  in  Applied  Mechanics  and  Engineering,  p.  269-289, 
August  1973. 

8.  Jones,  W.P.  and  Launder,  B.E.,  "The  Predictions  of  Laminarizat ion  with  a 
Two-Equation  Model  of  Turbulence,"  INT.  J.  Heat  Mass,  Transfer,  Volume  15, 
p.  68-87,  1972. 

9.  Launder,  B.E.  and  Spalding,  D.B.,  Lectures  in  Mathematical  Models  of 
Turbulence,  2nd  ed.,  Academic  Press,  1976. 

10.  Kays,  W.M.,  Convective  Heat  and  Mass  Transfer,  McGraw-Hill,  1966. 

11.  Boaz,  L.D.  and  Netzer,  D.W.,  An  Investigation  of  the  Internal  Ballistics 
of  Solid  Fuel  Ramjets,  Naval  Postgraduate  School  Report  Nunber  NPS-57Nt- 
7303U,  March  1973. 

12.  Mady,  C.J.,  Hickey,  P.J.  and  Netzer,  D.W.,  An  Investigation  of  the  Com¬ 
bustion  Behavior  of  Solid  Fuel  Ramjets,  Naval  Postgraduate  School  Report 
Number  N^-6^t^09^»  September  \§/1. 


35 


INITIAL  DISTRIBUTION  LIST 


1.  Library,  Code  0212 

Dean  of  Research,  Code  012 
Naval  Postgraduate  School 
Monterey,  CA  93940 

2.  Department  of  Aeronautics 
Code  67 

Naval  Postgraduate  School 
Monterey,  CA  93940 

Prof.  M.  F.  Platzer,  Chairman 
Prof.  D.  W.  Netzer 
LT.  C.  A.  Stevenson 

3.  Defense  Documentation  Center 
Attn :  DDC-TCA 

Cameron  Station,  Bldg.  5 
Alexandria,  VA  22 3 1 ^ 

4.  Naval  Air  Systems  Command 
Washington,  DC  20361 

AIR- 330 

5.  Naval  Weapons  Center 
China  Lake,  CA  93555 

Tech.  Library,  Code  753 
F.  Zarlingo,  Code  3246 
K.  Schadow,  Code  388 

6.  Chemical  Systems  Division 
United  Technologies 

P.0.  Box  358 
Sunnyvale,  CA  94088 
Tech.  Library 
R.  Dunlap 
A.  Holzman 
G.  Jensen 
P.  Willoughby 
P.  La Force 

7.  Chemical  Propulsion  Information  Agency 
APL-JHU 

Johns  Hopkins  Road 
Laurel,  MD  20810 

8.  AFAPL 

Wright-Patterson  AFB,  OH  45433 
R.  R.  Craig 
R.  D.  Stull 


No.  of  Copies 

2 

2 


1 

15 

2 

2 


2 


3 

3 

1 


2 


1 

I 


36 


